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We propose a mesoscopic setup which exhibits strong and manifestly non-Gaussian fluctuations 
of energy and temperature when suitably driven out of equilibrium. The setup consists of a normal 
metal island (N) coupled by tunnel junctions (I) to two superconducting leads (S), forming a SINIS 
structure, and is biased near the threshold voltage for quasiparticle tunneling, eV ~ 2A. The 
fluctuations can be measured by monitoring the time-dependent electric current through the system, 
which makes the setup suitable for the realization of feedback schemes which allow to stabilize the 
temperature to the desired value. 

PACS numbers: 74.50.+r,44.10.-|-i,72.70.-|-m 



In a (grand) canonical ensemble at temperature T, the 
total internal energy of the system fluctuates owing to the 
fluctuating energy and particle flows between the system 
and the thermal bath. It is a fundamental result of equi- 
librium statistical mechanics that the variance of energy 
is given by Var(£;) = kBC{T)T'^, C{T) being the (tem- 
perature dependent) heat capacity of the system [T]. For 
a system with a large number of degrees of freedom, the 
fluctuations are small and can be regarded as Gaussian. 

In many systems the internal relaxation rate is much 
faster than the rate of energy exchange with the envi- 
ronment. In this quasiequilibrium case, the probabil- 
ity distribution of the system is thermal with some ef- 
fective temperature T, which is unambiguously related 
to the instantaneous total energy of the system E via 
E/T = C(T). The fluctuations of energy and effective 
temperature are thus the same. In a driven system, the 
energy is not necessarily Boltzmann-distributed, so that 
the corresponding ensemble is non-Gibbsian. Moreover, 
the fluctuations of energy/temperature are not gener- 
ally Gaussian. However, in systems with a large num- 
ber of degrees of freedom these properties are usually 
non-accessible: The effective temperature is close to its 
average value determined from the heat balance, and its 
fluctuations are small and Gaussian. Recently, fluctua- 
tion statistics of effective temperature have been studied 
in non-interacting electron islands [5], and overheated 
single-electron transistors [3 2]. Typically, the non- 
Gaussian effects are noticeable only for large and there- 
fore exponentially improbable deviations from average 
values. 

In this Letter we demonstrate the feasibility of strong 
fluctuations of temperature and manifestly non-Gaussian 
distribution of these fluctuations in a mesoscopic system 
with a large number of degrees of freedom. The system is 
a SINIS structure, shown schematically in Fig. [ij where 
a normal metallic island is connected to two supercon- 
ducting leads via tunnel junctions, and biased close to 
the threshold for quasiparticle tunneling, eV « 2 A, A 
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FIG. 1. (color online) (top) SINIS structure biased at volt- 
age eV « 2A. The normal metallic island is connected via 
tunnel contacts to superconducting reservoirs. The effective 
temperature on the island, T, fluctuates due to fluctuations in 
the energy flows to the leads, (bottom) Energy diagram of the 
system for eV < 2A, showing the BCS density of states in the 
leads and the Fermi distribution on the island. The first or- 
der processes (single red arrow) are thermally activated. The 
second order processes (two blue arrows) occur at all tem- 
peratures and correspond to incoherent tunneling of Cooper 
pairs. 



being the energy gap in the superconductors. The cause 
of these fluctuations is the interplay of regular quasipar- 
ticle tunneling and two-electron Andreev tunneling. 

In this setup the temperature fluctuations can be easily 
and quickly monitored by measuring the electric current 
— no separate thermometers are necessary. This per- 
mits a practical realization of a feedback scheme where 
the fluctuations are coupled to certain control parameters 
of the system, so-called Maxwell demons [SHZ]- For ex- 
ample, when the system fluctuates to a low temperature. 
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the tunnel junctions could be "switched off," trapping the 
system at this temperature |8j. In one parameter regime 
discussed below, it is possible to reach an extremely low 
effective temperature this way. 

The basic physical mechanism responsible for the 
strong non-Gaussian fluctuations predicted here is the 
competition of quasiparticle and Andreev tunneling. Be- 
low and up to eV = 2A the quasiparticle tunneling pro- 
cesses cool the island, each tunneling event extracting on 
average an energy of fc^T ^ A. This is the well known 
cooling mechanism in SINIS structures [51-fTT]. Andreev 
tunneling, a process where a Cooper pair in the supercon- 
ductor is converted into two quasiparticles in the normal 
metal [12l [13] , deposits a relatively large amount of en- 
ergy, 2A, on the island. Energy and effective temperature 
are related hy E = Tr^fc^T^/ (6^), S A being the single- 
particle level spacing on the island, inversely proportional 
to its volume. A single Andreev event therefore heats 
the island to a temperature of at least ksTt = \/l2(5A/7r. 
Owing to the heat balance, the rate of Andreev processes 
r^i should be by a factor of fc^T/A smaller than the rate 
of quasiparticle tunneling. To this end, we may disregard 
the randomness in the quasiparticle flow and characterize 
the deterministic energy relaxation due to quasiparticles 
by a typical time r^. 

The regime of manifestly non-Gaussian fluctuations re- 
quires a low Andreev rate, Tj^Tj. < 1, implying |14| that 
all transmission eigenvalues Tp of the junctions should 
be small, J" = J2p'^p / J2p'^p ^ 1 being a crucial pa- 
rameter. This requirement is experimentally feasible: 
The Andreev events can be resolved in time [15], and 

~ 10~^ — 10~^ for aluminum junctions [TB]. More- 
over, the reservoirs should be kept at a low temperature 
and the island should be small such that it can be cooled 
down to temperatures of the order of Tt, bringing the 
average total energy of the island down to the order of A 
|17j . In addition, we need to avoid the Coulomb block- 
ade regime, so that the dimensionless conductance of the 



junctions should satisfy g = G/Gq > 1, Gq 



Finally, the heat exchange processes not involving elec- 
tron transfers should be small enough not to disturb the 
competition between Andreev and quasiparticle events. 

Under these conditions, each Andreev tunneling sub- 
stantially increases the temperature of the island. This 
increase is followed by the deterministic cooling at the 
timescale of t^, and the temperature remains low un- 
til the next Andreev event (see an example timeline in 
Fig. [2]) . The distribution of temperature is mainly deter- 
mined by the deterministic cooling and is strongly non- 
Gaussian. At an increased Andreev rate, V aTt ^ 1, 
the deterministic cooling is too slow to substantially de- 
crease the temperature between the Andreev events and 
the energy/temperature fluctuations become small and 
Gaussian. 

With these assumptions, the time evolution of the 
probability distribution function for the total energy on 
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FIG. 2. (color online) (top) Example timelines of energy (solid 
blue) and electric current (dashed red) for eV = 2A. Deter- 
ministic energy relaxation due to quasiparticle tunneling is 
characterized by a timescale Tr, relaxation time from E = 2 A 
to i5 = being Att for eV = 2A. Andreev tunneling hap- 
pens randomly at a rate V a- Current and energy are related 
by Eq. Q. (bottom) Voltage- Andreev rate "phase" diagram 
of the SINIS structure {u = {eV jl - /\)l(kBTt), 7 = V aTt)- 
Regions 1, 2, and 3 exhibit non-Gaussian fluctuations. 



the island, V(E)^ satisfies a Fokker-Plank equation 
dtV(E) - 



Be [H^^{E)V{E)^ 
Ta[V{E~2A)^V{E)]. (1) 



Here, the first term on the r.h.s. describes the deter- 
ministic relaxation due to quasiparticles, Hqp{E) being 
the energy flow of the quasiparticles, while the second 
term describes stochastic Andreev events with a rate F^. 
We assume a symmetric setup (TS] and a bias voltage of 
eV « 2 A ^ ksT. Under these assumptions, the energy 
flow reads \W\ 
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ikBTf/''A'/^F{iA - eV/2)/{kBT)), (2) 



F{x) - - 



rQ)Li3/2(-e--)+FQ)xLii/2(-e-^) 



F{x) is positive at x > 0, changing its sign at x ~ —0.72. 
Correspondingly, the quasiparticles always cool the is- 
land at eV < 2 A. At eV > 2 A, they fix the temperature 
of the island to ksT sa eV/2 - A [IT]. The Andreev rate 
is given bv [Tll^ 

F^ = gT-^ In [2A/ max(A - eV/2, keT)] . (3) 
47rn. 

It exhibits a weak logarithmic dependence on voltage and 
temperature that we disregard in the following. 
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Let us determine proper scales and corresponding di- 
mensionless variables in the parameter region of interest. 
The natural scale for the total energy is A, and we in- 
troduce a dimensionless energy, e = E/{2A). The island 
temperature in these units is given by T = Tt\fi^ and 
the natural scale for the bias voltage is correspondingly 
u — {eV /2 — / {kBTt) ■ In these units, the quasiparticle 
energy flow is 
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Tr 



(4) 



Dimensionless time is naturally expressed as t = t/Tr- 
The condition for manifestly non-Gaussian fluctuations 
discussed, TaTv = 7 < 1, can thus be expressed in 
terms of the tunnel parameter J^, ^ {ksTt/ A)^^^ ~ 
{5/ A)^/^. It is precisely the existence of this additional 
small parameter of the order of 6/ A which allows us to 
reach the regime of non-Gaussian fluctuations. 

The Fokker-Planck equation for the dimensionless 
variables reads 



drVie) 



(5) 



Its stationary solution depends on two dimensionless pa- 
rameters u and 7. There are six qualitatively different 
regimes, shown in Fig. [2] which we analyze briefly below. 

If the stationary distribution is Gaussian, we can ex- 
pand the difference V{e — 1) — 'P(e) in derivatives. In 
this case the average energy is determined from w(e) = 
^s/ipi-u/y^) = 7, and the variance can be estimated 
by Var(e) w (9c lnw(e))~^, Var(e) ^ 1 for a Gaussian 
distribution. There is a lower bound for the energy of 
the island where the cooling rate vanishes. It is given by 
the condition v{e) — with the result emin = (1.39u)^ for 
M > and emin = for w < 0. 

Let us first concentrate on the case |m| < 1, allowing 
us to approximate u/y^) « ^"(0) ~ 0.68. In the 
Gaussian case we have (e) — |Var(e) ~ (7/^(0))'*/'^. In 
terms of temperature (T) cx j'^^^A/kg, independent of (5. 
This is regime 6 in Fig.[2j The crossover to non-Gaussian 
behavior, regime 2 in Fig. [2j happens when (e) 1, i.e., 
7«F(0). 

In the case u 3> 1 the minimum energy is already larger 
than A. However, at sufficiently small 7 the distribution 
near this minimum energy is non-Gaussian. We can ap- 
proximate v{e) by expanding it near emin- Introducing a 
new variable e = e — emin we obtain a simpler equation 

Cu-^/^di [iV{i)] + 7 [P{f- - 1) - V{i)] = 0, (6) 



where C ~ 0.41. In the Gaussian case (regime 5 in Fig. [2]) 
we have (e) = Var(e) = u^^^^/C, with the crossover to 
non-Gaussian regime (regime 3 in Fig. [2]) taking place at 



7 w Cu~^/^, i.e., at a smaller value of 7 compared to 
the case |u| < 1. Qualitatively new behavior takes place 
when (e) ~ emin, that is, when 7 > u^^'^. In this case 
we can approximate F^—u/y/e) ~ F{Q), and recover the 
same behavior as in the case |w| ^ 1 (regime 6). 

Finally, in the case |u| ^ 1, u < 0, the quasiparticle 
rate is exponentially suppressed at low energies, i/e < 
Let us assume that the distribution is concentrated near 
some value of energy e^ <C u^. Using the approximation 
F{x) ^/iixe^^ for x 00, we have near e e* 



V^|«|ei/^e-l"l/^9, 



+ ^[V{e-l)-V{e)]=Q, (7) 

where e is the deviation from e* . In the Gaussian case, we 
can estimate the variance as Var(e) ~ 2e*^^/|u|. Com- 
paring this with unity, we conclude that the Gaussian 
distribution is realized if e* > \u/2\'^/^ <^ |u, i.e., if 
7 > |u|^/6exp(-2i/3|y|2/3) (regime 4 in Fig. 2|. For 

smaller 7 the distribution is non-Gaussian (regime 1 in 
Fig. [2]). Upon further increase of 7, e* grows, allowing us 
to again approximate F{—u/y/e) F{0). It reaches 
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, taking us back to regime 6. This is similar 



to the situation with m 3> 1. 

The probability distributions for energy, calculated nu- 
merically from Eq. ^ , are shown in Fig. |3] We choose 
three different values of u corresponding to the three non- 
Gaussian regimes, and vary 7 to observe the crossover to 
Gaussian behavior. 

The topmost figure shows the probabilities for u = 0. 
The distributions for the two smallest values of 7 are 
strongly non-Gaussian, belonging to regime 2 in Fig. [2] 
For larger 7 the average energy grows, being proportional 
to 7"^/^. The inset in the figure shows some realizations 
of energy timelines corresponding to the distributions. 
For small 7 one can clearly see the difference between 
Andreev events, which are Poisson distributed and in- 
crease the energy by 2A, and quasiparticle processes, 
which steadily cool down the island. For large 7 the 
energy-increasing and energy-decreasing processes look 
the same, characteristic of a Gaussian distribution. As 
seen from the timeline, the energy on the island can reach 
zero and stay there for a finite time, implying that the 
temperature goes to absolute zero. In reality, the tem- 
perature would be about S/kg <C Tt when only a few 
particles are excited on the island, a situation outside 
the scope of the present model. 

The middle figure corresponds to u = 1. Since u > 0, 
the minimum energy is finite, Cmin ~ 1.93. The distribu- 
tion for the smallest value of 7 is clearly non-Gaussian, 
belonging to regime 3 in Fig. [2] For lar ger 7 the average 
energy grows linearly. 

The last figure shows the probabilities for u — —2. 
Again, the distribution for the smallest value of 7 is non- 
Gaussian, belonging to regime 1 in Fig. [2] Upon an in- 
crease of 7 the main body of the distribution transforms 
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FIG. 3. (color online) Probability distributions for the inter- 
nal energy on the island for (top to bottom) u — 0, u = 1, and 
It — —2, and some values of 7. Insets show realizations of en- 
ergy timelines for each value of 7. The bins of the histograms 
have a width of 0.1. 



into a Gaussian one and the average energy grows. The 
tails of the distribution deviate from Gaussian even for 
the largest 7 shown. 

Since the electric current through the SINIS structure 
is sensitive to the instantaneous temperature on the is- 
land, the statistics of internal energy and temperature 
can be gathered by monitoring the instantaneous electric 
current. The necessary relation is given by 



Tr KbTi \ 2 



which can be numerically inverted. The time scale of 
the fluctuations, Tr, is of the order of 10 ns, and the 
corresponding current scale, eA/{TrkBTt), of the order 
of 1 nA for S = lO'^ K x fc^, A = 1 K x fc^, = 10 mK 
and g = 5, corresponding to a copper island of volume 
V = 0.01 (^m)^ connected to aluminum leads. 

In summary, we suggest a prototypical setup utilizing a 
voltage biased SINIS tunnel structure for the detection of 
energy and temperature fluctuations out of equilibrium. 
We have identified parameter regimes which exhibit fluc- 
tuations that are both strong and non-Gaussian, and we 
suggest that the measurement of instantaneous electric 
current could be used to gather the statistics of these 
fluctuating quantities. We have also described a possibil- 
ity to realize Maxwell-demon-like feedback schemes. 
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